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SUMMARY 

A new generalization of Dawson's integral function based on the link between 
a Riccati nonlinear differential equation and a second-order ordinary differ- 
ential equation is reported. The MacLaurin expansion of this generalized 
function is built up and to this end an explicit formula for a generic cofactor 
of a triangular matrix is deduced. 

1. INTRODUCTION 

The study of homogeneous linear second-order ordinary differential equa- 
tions (HLIIODE) is strictly related to that of Riccati equations. It is possible 
to express the general integral of an arbitrary HLIIODE in terms of a par- 
ticular solution of the associated Riccati equation. In the first part of this 
paper we recall this procedure in order to write down explicitly the most gen- 
eral solution of a HLIIODE vanishing at a conveniently chosen fixed point. 
The structure of such a general solution suggests a very simple way to extend 
Dawson's integral function (DIF), appearing in various physical contexts such 
as spectroscopy, electrical oscillations, heat conduction, astrophysics and so 
on as well as in applied mathematics [T]-|10j. Quite recently the DIF has 
emerged in the treatment of the dynamics of the so-called generalized spin 
star system where the reduced dynamics of a system composed of two cen- 
tral qubits is investigated in the limit of an infinite number of environmental 
spins [TTJ. The DIF deserves, in addition, attention on the mathematical side 
since it is related to other special functions and even because it possesses use- 
ful formal properties exploitable for its numerical evaluations [121 US]- The 
main scope of this paper is to propose a new class of transcendental functions 
which may be viewed as generalizations of the DIF. We report in detail the 
construction of the MacLaurin expansion of a generic element of this class 
by solving a linear system of infinitely many equations in infinitely many 
unknowns. We reach this goal evaluating all the cofactors of a triangular 
matrix of arbitrary finite dimension having all its diagonal elements equal to 
one. 
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2. RICCATI GENERAL SOLUTION OF A HLIIODE 



Consider the following homogenous second-order linear differential equa- 
tion with variable coefficients 

y"{x) + b l {x)y'{x) + b 2 {x)y{x) = (I) 

where b\ (x) and b 2 {%) are C 2 functions on (—a, a) e R, a > 0. Such a linear 
equation is invariant under dilatations y 1— > A?/ the infinitesimal generator of 
which is 

Lie theory of symmetry of differential equations takes into account such sym- 
metry and, as indicated in [H], the recipe is to look for a new coordinate u 
such that X becomes 

The coordinate X is such that Xu = 1, i.e. y|| = 1, from which we obtain 

y(x) = e u ^ (up to multiplication by a constant). 
Since in terms of the new coordinate we get 

y'{x) = u'(x)y(x), y"( x ) = u "{ x ) + y{ x ) (4) 

the transformed differential equation may be cast in the following form: 

u"{x) + (V(x)) 2 + bi{x)u\x) + b 2 {x) = (5) 
which takes the form of a Riccati equation 

z'(x) + z 2 (x) + b 1 {x)z{x) + b 2 (x) = (6) 

if we put 

z(x) = u'(x). (7) 

Let z(x) be a particular solution of eq.©. Then on integrating eq.([7]) 
one immediately gets the particular solution of eq.(pQ) 

yi (x) = e%*® dt (8) 

such that yi(0) = 1 and y[(0) = z(0). 
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We recall that, once a solution is known, we can introduce the change of 
variable y(x) = y\(x)v{x) and then, taking into account that 

y'(x) = y[(x)v(x)+yi(x)v' (x) y"{x) = y'[ (x)v(x)+2y[(x)v' (x)+yi(x)v" (x) 

(9) 

and the fact that yi(x) is a solution of eq.([T]), eq.([T]) becomes 



v"(x) = a(x)v'(x) 



(10) 



where we have defined 



A solution of eq.l fTOi) is 



a{x) = —^2^(2;) + b\(x) 



v(x) 



exp 



J a{t')dt'\dt 



from which we obtain the following solution of eq.flT]): 







/ exp 


[ a(t')dt' 


/o 


Jo 



dt 



satisfying the initial conditions ?/2(0) = and y' 2 (0) = 1. 
In view of the fact that 



W(0) = 



2/i(0) y 2 (0) 

y'M 2/2(0) 



1 

z(0) 1 



1^0, 



we immediately write down the following general solution of eq.(pQ): 



(12) 



(13) 



(14) 



y[x) = exp 



z{t)dt 



Uo 



C 1 + C 2 / exp 



2z{t') + b 1 {t'))dt' 



dt 



_ (15) 

Requiring y(0) = yields C\ = so that the most general solution y(x) of 
eq.Q vanishing at x = may be cast in the form 



y(x) = C exp 



z{t)dt 



exp 



2z(t') + b l {t'))dt' 



dt (16) 



which satisfies the further Cauchy condition y'(0) = C. 
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3. EXTENDING DAWSON'S INTEGRAL FUNCTION 



The form of y(x) provides a very favorable starting point to construct 
generalizations of Dawson's integral odd function 

PX 

F(2,x) = exp(-x 2 ) / exp(t 2 )dt = -F(2, -x) (17) 
Jo 

where i6l. In the literature there have appeared extensions like 

F(p,x) = exp(-x p ) [ exp(t p )dt (18) 



with p = 2,3,4,... and x G R, in connections with interesting problems 
both in applied Physics and in Mathematics [T5J, [16] . Our scope is to exploit 
the structure of y(x) to propose new generalizations of the original Dawson's 
integral function. The first natural extension suggested by equations (ITT]) and 
(115)1 is to look for bi(x) in equation (11) such that 

z(x) = 2z(x) + bi(x) z(x) = -bi(x) (19) 

This condition establishes the following link between bi(x) and b 2 {x) through- 
out the associated Riccati equation (Jf"} 

- b[{x) = -b\{x) + b\(x) - b 2 {x) b 2 (x) = b[(x) (20) 

so that we may claim that the unique solution of the Cauchy problem 



y"(x) + b(x)y'(x) + b'{x)y{x) = y"{x) + [b(x)y(x)J = ^ 
2/(0) = y'(0) = l 

may be written down as 

y(x) = exp[-B(x)] / exp[B(t)]dt (22) 
Jo 

where 

B{x) = / b(t)dt. (23) 



When b(x) = 2x or b(x) = px p_1 eq. (l2*2l gives back the function F(2, x) and 
F(j>, x) respectively. We denote the function defined by eq. (122*1) . by D b (x) and 
call it the generalized Dawson's integral function associated to the function 
6(1). 
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There is another way of extending Dawson's integral function easily sug- 
gested by equation ffTBj) . If we in fact stipulate that 



z( x ) = -XpxP' 1 pGN+AeR-{0} 
2z(x) + bi(x) = -fisx 3 - 1 s G N + , /i£K - {0} 



(24) 



then b\{x) = 2Xpx p 1 — fisx s 1 , and consequently 

b 2 (x) = X 2 p 2 x 2p - 2 - Xfipsx p+S - 2 + Xp{p - l)x p - 2 . (25) 
Thus we may claim that the function 

F[(X,p), (n,s);x\ = exp(-Xx p ) / exp(fit s )dt (26) 

Jo 

with p, s G N + , A,/iGR- {0}, iGlis the unique solution of the following 
Cauchy problem 

y"{x) + f^Xpx 1 "- 1 - fisx s - 1 ^jy'{x) + 

+ (x 2 p 2 x 2p - 2 - Xiispx s+p - 2 + Xp{p- l)x p - 2 ^y(x) = • ( 27 ) 
y(0) = 1/(0) = 1 



F[(l,p), (l,p); x] is of course coincident with F(p, x) and is included in eq.(!22 
in correspondence with 6(x) = px p ~ l . 



4. MACLAURIN EXPANSION OF D fe (x) 



We concentrate on D^x) posing the following question: if b(x) may be 
expanded in a MacLaurin series with the radius of convergence R, that is 

n=0 

&( n )(:r) being the n-th derivative of 6(x), is it possible to find the MacLaurin 
expansion of D^(x) in terms of the class of coefficients |£/ n )(0), n G Nj ? 
The question is well posed, since the assumption on b(x) guarantees that 
Db(x) too may be expanded in MacLaurin series with the same radius of 
convergence R. Then our problem is to build up the explicit expression of 



5 



D b k \o) as a function of {6^(0), n G N}. To this end we begin by observing 
that for any b(x) 

D' b (x) = -D b (x)b(x) + 1, D'ltx) = - [D b (x)b(x)] ' (29) 

D b (0) = 0, D' b (0) = 1 (30) 

Dl™\x) = -[D b (x)b(x)f +1) = -f:( k+1 )Di n \x)b^\x) (31) 

n=0 ^ ' 



Df +1) (0) = -kb^iO) - ( J£>i n) (0)6 (fc - n) (0)- 



n=2 



(32) 



To find "0(0), with A; G N + , we must solve the linear system of k equations 
in the k unknowns j-D^(O), D^(0), . . . , D^ k+1 ^ (0) j . Its incomplete matrix 



is 



( 1 

©6(0) 1 

©6«(0) ©6(0) 1 

gW(0) ©6W(0) ©6(0) 1 

©6(3)(0) ©6^(0) ©6«(0) ©6(0) 1 

V ©6^(0) ©M fc - 3 )(0) ©6( fc ~ 4 )(0) 



\ 



©ko) iy 

(33) 

and has lower triangular form with all its diagonal elements equal to one. 
The vector Bk of the constants appearing in the linear system (1521) is 



( -6(0) \ 
-Q6«(0) 

V -©6( fcLi )(o) y 



(34) 



It is convenient to introduce the following notation for the elements of of 
Ah and 6j of 



">,3 



i < j 

1 z = j 

oii)6 (< " i_i) (o) z>j 



^-^(O) = -26 (i - 1} (0) 



(35) 
(36) 



provided that 1 < i,j < k. Since det = 1 for any k G N + , Cramer's 
theorem yields 

Df +1) (0) = detP fc (37) 



where 



Ph 



1 

g)6(0) 1 

©&W(0) Q6(0) 1 



-6(0) \ 
-26«(0) 
-36( 2 )(0) 



(38 



V g)6 (fc " 2) (0) g)M fc - 3 )(0) ©6(0) -fcft^-^CO) / 



The determinant associated to the matrix borders the determinant of 
A^_i through the addition of the fc-th row and the fc-th column. The value 
of det Pfc may thus be evaluated by means of the well known Cauchy formula 
[T7] according to which 



det B = det 



Oil ■■■ 

Cnl ■■■ 

V Pi ft 



«1 \ 

a- 2 



fin CL J 



a det C — a r /? s A riS (39) 



r,s=l 



where A r >s is the cofactor of the corresponding element c rjS of the bordered 
matrix 

/ Cn C12 ... Ci n \ 

c 2 i 



C 



(40) 



/ 



\ Cjii c nn 

In our case B = P^ and C = A*._i, so that this formula assumes the form 



/j(* +1 )(0) = detP fc = -fc6(*-^(0) + ^ if . U^- 1 )(0)^- 1 )(0)A i ,(41) 

ij=i v + / 

Ajj being the cofactor of the correspondent element ay (1 < < A; — l)of 
the matrix Ak-i- It is easy to convince oneself that Ajj = for any i > j, 
since the elimination of the i-th row and j-th column with i > j gives a 
triangular matrix with one diagonal element equal to 0. It is in addition 
evident that A^ = 1 for any i. Thus we have only to evaluate the cofactors 
C(i,n) = Ai^ +n with 1 < i < k — 2 and I < n < k — i — 1 finding the explicit 
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expression of such cofactors for a generic lower triangular matrix C = (c rjS ) 
of order k — 1 with all its diagonal elements equal to 1. We indeed claim that 

n 

C(i, n) ^ ( 1) ^ Cj+ n j_)_p :L Cj_|_p 1) j-|_p 2 . . . Cj_)_p s _ 2i j_|_p s _ 1 Cj_|_p s _ 1) j_|_p s (42) 

8=1 {p}? 

where denotes En>p 1 >p 2 >-> Ps . 1 > Ps =o tliat is summation over all the 

(nlD possible sets of s indices {pi,P2, • • • ,Ps}, taken in strict decreasing order 
and such that < Pi < n, i = 1, 2, . . . , s — 1 and p s = 0. We observe that 
the terms of eq.( H2i) corresponding to a fixed value of s are products of s 
elements of the matrix A with a plus (minus) sign if s is even (odd). Then, 
the only term corresponding to s — 1 is —c i+n ^. 

Let us prove eq. (1421) by induction on n for an arbitrary i = 1, 2, . . . , k — 2. 
Firstly we show that 

C(i, 1) = -Ci + i,i. (43) 

To obtain C(i, 1) = A^+i we have to calculate and multiply by —1 the 
determinant of the matrix obtained from Aj s _i eliminating the z-th row and 
the (i + l)-th column. It is easy to observe that this operation leaves a 
triangular matrix having q+i^ as a diagonal element, and all other diagonal 
elements equal to 1. As a consequence, we have A iii+ i = — Cj+i,i, as stated 
before. Let us now suppose that, for all m < n we have 

m 

C(J,,Ul) = ^ { 1) ^ ^ Gi+m,i+piGi+pi,i+p2 ■ ■ ■ c i+p a -2,i+p s -i C i+Ps-i,i+Ps • (44) 

8=1 {pyn 

Using Laplace's second theorem with the i-th row for the cofactors and the 
(i + n + l)-th row for the elements we obtain 

k-l 

q+n+ijAjj = 0. (45) 

3=1 

Using the fact that C(i, 0) = 1, cjj = for j > % whereas = 1 we obtain 

i+n 

Cj+n+i.i + c *+™+iJ C '(^i ~~ + C '(^' " + 1) = 0. (46) 
j=i+i 

Using eq. (I44p and posing j — i = m we have 



C(i,Tl + 1) — — Cj+n+i^ + ^ ] ^^( — 1) S+1 ^ ] Cj+n+i^j+mQ 



j+m,i+pi • • • c i+p 3 -i,i+ps 



m=l s=l {p}j 



Observing that for a fixed value of s the second term contains products of 
s + 1 elements qj with the correct sign factor, it is easy to convince oneself 
that we may write 

n+1 

C(i, n + 1) — —C i+n+ i :i + ^^( _ I)* Ci+n+l,i+piCi+pi,i+p 2 • • • Q+p 3 _i,i+p s (48) 

s=2 {p} (n+i) 

that is 

n+1 

C(i, n + 1) = ^^(— l) s c i+n+ i )i+pl Cj + p 1)i+ p 2 . . . c i+ p s _ 1>i+ p s (49) 

3=1 M«" +1) 

which concludes the demonstration eq.( j4^|) . 

As far as we know the explicit expression of all the cofactors of a triangular 
determinant of arbitrary finite order having all its diagonal elements equal 
to one has not previously appeared in literature. 
Summing up we have proved that, for any 1 < i, j < k — 1, 

(C(i,j-i) i<j 
1 1 = 3 . (50) 

i>j 

This concludes the derivation of the MacLaurin expansion of Db(x) for any 
b(x) as well expandable in MacLaurin series. We guess that other general 
properties may be proved for such class of extended Dawson's integral func- 
tions also thanks to the explicit knowledge of their MacLaurin expansions. 
Moreover special cases obtained in correspondence to particular choices of the 
function b(x) may find applications both in Physics and in Applied Mathe- 
matics. 

5. CONCLUDING REMARKS 

The structure of the general integral of a linear homogeneous second-order 
ordinary differential equation in the form given by eq. (fT5"]) suggests a very 
simple way to generalize Dawson's integral function. Such a generalization 
introduces indeed a class of new functions all possessing the same Dawson - 
like structure. 

In order to explore some properties of this class without choosing the form 
of b(x), we have coped with the interesting question of finding the MacLaurin 
expansion of D\,{x) in terms of the MacLaurin coefficients of the same b(x). 
It is worth noting that the resolution of this problem has lead us to derive 
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the explicit expression of all the cofactors of a lower triangular matrix which 
diagonal entries are all equal. We underline that such a derivation may 
be successfully generalized and applied to a generic triangular matrix. In 
conclusion we feel that the ideas and methods reported in this paper may 
be of some help to propose new generalizations of other special functions of 
interest in physics, chemistry and applied mathematics. 
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